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We consider a model of neutrino oscillations with three additional sterile neutrinos and a gauged 
B — L interaction. We find allowed values of the model parameters which can reconcile the results 
of the evidence for conversion seen at the LSND neutrino oscillation experiment with 

the null results of the MiniBooNE experiment. A portion of the low energy excess of Ve. events 
seen at MiniBooNE can arise naturally, and we make a quantitative prediction for the forthcoming 
anti-neutrino oscillation results at MiniBooNE. 



I. INTRODUCTION 



The standard three neutrino oscillation model is experimentally established by long baseline neutrino experiments 
[TJIHISlllllSllSllZllSlini HO] , which are sensitive to differences in neutrino mass squared which are much smaller than 0.1 
eV^. The neutrino flavor eigenstates e, /z, r are related to the mass eigenstates via a 3-by-3 unitary matrix, neglecting 
a CP-violating phase 8: 

Cl3 Si3 \ / Ci2 Sl2 
10 -Si2 Ci2 I I I , (1) 
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where Cij — cos Qij , = sin Bij . The three angles Qij have been measured to good accuracy: 

tan^ 6*12 = 0.45 ± 0.05 

sin2 20i3-Otr' (2) 
sin2 2023 = l.Oi;^.i 

The neutrino mass squared differences have also been determined: 

Am?2 = (8.0 ±0.3) • 10-5 eV^ 
lAm^gl = (2.5 ± 0.2) • 10^^ eV^ 
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While many experiments have confirmed this picture, a pair of short baseline experiments suggest new physics may 
be contributing to neutrino oscillations. The LSND experiment found evidence at the 3cr level for v^^ oscillations 
[TTJ [121 US]- This result does not fit into the standard oscillation model, as the LSND experiment was sensitive to 
oscillations between neutrinos with a mass squared difference greater than 0.1 eV^. 

In an attempt to confirm or rule out the LSND results, the MiniBooNE experiment searched for ly^ oscillations 

over the same Am^ range as LSND, but with higher energy and longer baseline. The initial results of the experiment 
are inconsistent with a standard oscillation interpretation of LSND, but saw an excess of electron candidates in the 
low energy spectrum of their data [14) . The low energy excess occurs at energies below 475 MeV, in a region not 
included in the oscillation analysis, and where the uncertainty in the background is larger than in the analysis region. 
If the excess at MiniBooNE is a physics signal rather than some nonunderstood background, then new physics is 
needed, which might also be able to explain the results of LSND. Several attempts have been made to reconcile 
the LSND and MiniBooNE results with each other and with other oscillation experiments. This appears to require 
new ingredients besides simple neutrino oscillations, even if more sterile neutrinos and CP violation are introduced 

[iniiiiiiiniiiTiiiHiiisiEo]. 

To accommodate the results of the MiniBooNE and LSND experiments, we introduce a model for neutrino oscil- 
lations that can be compatible with the long baseline results, the various null short baseline results, and the LSND 
and MiniBooNE anomalies. As in the conventional seesaw model [211 [22l |23l [24l |25l [26] , neutrino masses are allowed 
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via the introduction of 3 sterile neutrinos. However, the sterile neutrinos have small masses in the eV range. (Such 
"miniseesaw" models have been considered before, for a variety of reasons |27, 28, 29, 30, [311 131].) Vacuum mixing 
angles for oscillations between the active and sterile sectors are all fairly small. A large MSW-type potential for 
neutrinos in matter |33} [34], arising from a, B — L gauge interaction, gives strong energy dependence in the short 
baseline small mixing angle oscillations. Three sterile neutrinos are an integral part of the model, as they are required 
for gauge anomaly cancellation and to allow for neutrino masses. Long baseline data indicates that new interactions 
for the active neutrinos should either be flavor diagonal or weak |35l 1361 137] . However because the active neutrinos 
all have the same B — L charge, in the absence of mixing with sterile neutrinos the B — L potential has no effect on 
the oscillations among active neutrinos. While the B — L potential does give a significant matter effect for the mixing 
with sterile neutrinos, when the active-sterile mixing is small, the effects of the new interaction on long baseline, large 
mixing angle oscillations among active neutrinos can be negligible. 

We note that in order to get a significant matter effect at short baseline, we will need the B — L boson to be light. 
We will argue that there is a phenomenologically allowed window for a weakly coupled U{\)b~-l gauge interaction 
with a light (tens of keV) vector boson with a coupling constant less than ~ 10~^. While the effects of such a boson 
are small in the precision electroweak and QED experiments, there can be a significant effect on neutrino oscillations. 

Several features of the LSND and MiniBooNE experimental results arise naturally in the model. Short baseline 
active flavor oscillations at higher neutrino energy are suppressed by the MSW-type potential in matter, and anti- 
neutrinos tend to have higher oscillation probabilities than neutrinos. The results of MiniBooNE relative to LSND 
can be qualitatively explained by these features - MiniBooNE searched for oscillations at higher energy and with 
neutrinos, while the strongest evidence for oscillation found at the lower energy LSND experiment was with anti- 
neutrinos. Further, part of the low energy electron neutrino excess at MiniBooNE can be explained by the strong 
energy dependence of oscillation parameters in this model, although it is difficult to fit the entire excess. 

After discussing the model and deriving the formula for oscillations, we proceed to fit the results of LSND, Mini- 
BooNE, and short baseline neutrino disappearance experiments numerically. We find a substantial portion of param- 
eter space capable of fitting all experiments simultaneously, and give a prediction for the outcome of the MiniBooNE 
anti-neutrino oscillation experiment. 



We give a brief summary of the two experiments. LSND observed a 3.2 cr excess of oscillation events using 

jX^ decays at rest, which produced a muon anti- neutrino energy spectrum up to the Michel endpoint for muon decay. 
A secondary analysis looked for oscillations in from 7r+ decays in flight, which yield neutrino energies in the 

range 60 to 200 MeV. This analysis initially observed a 2.6 a excess, but a later analysis showed a reduced signal of 
0.6 a. The observed anti-neutrino oscillation signal is 



The LSND detector was centered 30 m from the neutrino source, meaning the anti-neutrino oscillation experiment 
was sensitive to a neutrino mass squared difference 



MiniBooNE was designed to be sensitive to the same range in Am^, but is characteristically quite different from 
LSND. MiniBooNE searched for — > v^. oscillations using muon neutrinos with energies between 200 MeV and 3 
GeV, and the detector was situated 541 m from the neutrino source. 

In their recently released results, MiniBooNE found no evidence for neutrino oscillations in the energy range 475 - 
3000 MeV, but an unexplained low energy excess of Vf. in the energy range 200 MeV to 475 MeV. The fit of the excess 
to conventional oscillations is poor, and the excess corresponds to percent level oscillations in this energy regime. The 
result of an analysis of oscillations at MiniBooNE is expected in the coming months, and could prove to be 

a strong test of the anomalies seen at both LSND and MiniBooNE. 

III. A RENORMALIZABLE MODEL OF NEUTRINO MASSES WITH A LIGHT B-L GAUGE BOSON 

We motivate the model by considering the phenomenology of a flavor universal MSW-type matter effect. We take 
the results of LSND and MiniBooNE at face value: they found different oscillation results, and the key differences 
between the two experiments are the energy scales and the use of neutrinos at MiniBooNE and anti-neutrinos at LSND. 
We thus need an energy dependent effect and/or one which distinguishes between neutrinos and anti-neutrinos. Vector 
interactions with matter do both. We take this as evidence for a new gauge interaction, whose potential provides a 
strong energy dependence in short baseline oscillations. In this model we utilize a weakly coupled, flavor universal 
U{\)b~l gauge interaction to generate such a potential. The effects of such a flavor universal interaction are different 



II. THE LSND AND MINIBOONE EXPERIMENTS 




De) = (0.264 ±0.081)% 




(4) 
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from those of flavor dependent potentials which have been considered before [331 IMJ [3SJ ISZl [3HJ ISHJ SOI SIJ 112] , in 

particular, the effects on long baseline active flavor conversion are small. In order to obtain a sufficiently large effect 
to impact short baseline oscillations, we require a new light gauge boson with mass in the several keV region. 

We extend the standard model to include the B — L gauge interaction by making all derivatives covariant under 
a, B — L gauge interaction as well as under the standard model gauge group, and adding the conventional kinetic 
term for the B — L gauge boson. We refer to this gauge boson as the paraphoton. Gauge anomaly cancellation is 
ensured by the addition of 3 neutrinos which are sterile under the standard model interactions and have opposite 
B — L charge to the usual neutrinos. As discussed further in §VI A| the paraphoton acquires a small mass from a 
condensate of a scalar field s. A Yukawa coupling to the condensate will be responsible for the Majorana masses of 
the sterile neutrinos. 

Because neutrino energies far exceed their mass in oscillation experiments, we can expand the effective Hamiltonian 
for neutrino propagation: 

n^E-l + ^^+V (5) 

where E is the neutrino energy, M is the mass matrix, and V is the potential matrix. V receives contributions both 
from the standard model weak interactions and from the B — L interaction. 

Neutrino masses can be Dirac {B — L conserving) or Majorana {B — L violating). Both types of mass terms are 
required to give active-sterile mixing. Using a notation where all fermion fields are 2 component left handed Weyl 
fields, and suppressing the Lorentz indices, the allowed Dirac mass terms are 

niijViNj + h.c. , (6) 

where we have labeled the sterile neutrinos Ni, i — 1,2,3. These mass terms can arise from a small Yukawa coupling 
to the usual Higgs boson. Majorana masses for the sterile neutrinos arise from a coupling 

y^jsN.Nj + h.c. (7) 

where yij are dimensionless couplings, and s is a scalar field with B — L charge -2 and no standard model gauge 
interactions. Then the sterile neutrino Majorana mass matrix will have entries 

= y^{s) (8) 

and their mass terms are 

M,jNiNj+ h.c. (9) 

The fiavor basis for sterile neutrinos is irrelevant, so we take the sterile Majorana mass matrix A^s to be a 3 by 3 
diagonal matrix, with eigenvalues M^, for i — 1^2, 3. Therefore, the mass matrix for all 6 neutrinos is 

where AId is the 3 by 3 Dirac mass matrix. Atmospheric and solar oscillation experiments provide good constraints 
on large mixing angle long baseline active fiavor mixing. As we will see, in matter at high energy, the active neutrinos 
nearly decouple from the sterile neutrinos, and mainly oscillate between active fiavors with parameters governed by 
the Dirac mass matrix. Therefore, we will take the mass squared differences in the Dirac mass matrix to be small. 
This makes the Dirac masses almost degenerate, with an overall scale m, which we call the Dirac mass from now on. 
Then the Dirac mass matrix can be written 

Md ~ mu (11) 

where u is a unitary 3-by-3 matrix relating the mass and fiavor bases. Long baseline tests are sensitive to small 
non-degeneracy in the masses, while short baseline tests are insensitive. Parameterizing Md = mu + 5m, with the 
entries in 5m much smaller than to, we note than 5m may be diagonalized in a different basis from mu, and that small 
5m has a significant effect on long baseline oscillations but a negligible effect at short baseline. The mixing angles 
relevant for long baseline depend both on the angles in u and on 5m, while the parameters relevant for short baseline 
disappearance and active fiavor conversion depend on u, m and the three Majorana masses Mi. Because long and 
short baseline oscillations depend on independent parameters, we do not consider the constraints from long baseline 
oscillations in our analysis, and neglect 5m. 

We do not wish to do detailed model building in this paper, but note that a nonabelian fiavor symmetry and a 
suitable choice of fiavor breaking spurious could explain the near degeneracy in the Dirac mass matrix. 



Note that, as discussed further in [VIA the expectation value of s may be affected by the environment 



l46l l47l l48l l49l ISOj [5TJ [52l |53], and so the paraphoton mass and neutrino oscillation parameters could be quite 



different in extreme environments such as stars, supernovae, and the early universe. For this reason the constraints 
from astrophysics and cosmology [SHllMlfSS, 56, 57, 58, 59i [301 ISD EU [S3] are compHcated and depend on different 
parameters than the neutrino oscillations. We assume that the neutrinos masses are constant in our analysis of 
terrestrial experiments, although it could be interesting to consider environmental dependence of the neutrino masses 

miiHiiMi. 

We can now write the mass and potential matrices in the flavor basis in terms of 3-by-3 blocks: 

We rotate from the flavor basis to a basis where the Dirac mass matrix mu in Ad is diagonal, via the unitary matrix 
U, which leaves V unchanged: 

We call this the prime basis. In this basis the effective Hamiltonian for neutrino evolution is, 

^ ^ 2E 2E\ mMs 2VE-t + Ml)- ^ > 

Note that we are neglecting the usual MSW terms from the weak interactions, as these have negligible effect at short 
baseline. The Hamiltonian in the prime basis is 2-by-2 block diagonal, with each block having the form 

-w' - p j_ ^ ( ^ ^ (^ K\ 

- ^ ^ 2E ~ ^ ^ 2E [mM, WE + Mf ) ' ^^^> 

If the sign of V is positive, then provided that m <C Mi, we have a mini seesaw with small active-sterile mixing. 
In this case the lighter neutrinos are almost entirely active, while the heavier are almost all sterile. Furthermore, the 
probability of active neutrino disappearance into sterile is suppressed by the ratio m? /Mf . As the energy of the ex- 
periment increases, the active-sterile mixing is suppressed further. This high energy limit is relevant for MiniBooNE, 
which has an energy that is an order of magnitude larger than LSND, and can therefore have highly suppressed oscilla- 
tions. Note that the small mass differences in the mostly active neutrinos contribute essentially nothing to oscillations 
at MiniBooNE and LSND, just as the standard oscillation model predicts no oscillations at these experiments, but 
will be important when considering long baseline oscillation experiments. 

If the sign of V is negative, as is the case for anti-neutrinos, then there can be a range of resonance energies for 
which AVE ~ Mf , giving large mixing between the active and sterile sectors. This enhances oscillations between 
active flavors, and may account for oscillations seen in anti-neutrinos at LSND. Oscillations are suppressed when 
|4Fi? -I- Mf\ ^ mMi, so the effects of sterile neutrinos on anti-neutrino oscillation experiments should be small at 
high energy, but larger than on neutrino oscillations. 



(13) 

V = V 



IV. OSCILLATION PROBABILITIES 



Here we flnd the formula for oscillation between active flavors. We pay special attention to those experiments 
where the potential could be relevant, and discuss the implications of solar and atmospheric neutrino experiments. 
We outline the constraints on the B-L gauge interaction parameters from electroweak measurements and discuss other 
possible constraints. Finally, we fit to the results of LSND and MiniBooNE, and make a prediction for anti-neutrino 
oscillations at MiniBooNE. 

To find the oscillation probability, we diagonalize the Hamiltonian and find the transformation between the fiavor 
basis and diagonal basis. In the prime basis, we utilize the block diagonal nature of the Hamiltonian. Each 2-by-2 



block in (14) is diagonalized via the unitary angle 9i, where 



2mM, 

tan 29^ = (16) 

These angles define the unitary matrix needed to diagonalize the Hamiltonian in the prime basis: 

Ud =u''^n'u' 

to2 1 (17) 

= S + ^ - ^ + ^diag(Ar, A^, A3-, A+, A+, A3+) 
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AVE + Mf 



2V ' ^Am^Mf + {AVE + Mf)^ , 
Between the flavor basis and the diagonal basis, we have made the transformation 



(20) 



The oscillation probability for neutrinos of a specific energy E over a baseline L is given by the standard oscillation 
formula 1651. 



Pa^b{L, E) = ^ ^Kab,af3 COS (^AapL - a.rg(Kab,af3 

a, 13 



where 



Kab.ap — I^d" ^'d* ^[f* ^IJ 



(21) 
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For simplicity in the numerical analysis we will use the values ^13 = and 623 = 7r/4, the angles suggested by long 
baseline oscillation experiments, although in principle, the angles need not be related. The probability of oscillation 
between electron and muon neutrinos is 
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Before moving on to numerical analysis, we discuss constraints on the parameters governing the oscillation probability 
from experiments. Long baseline oscillation tests help constrain the parameters affecting the oscillation probability. 
Experiments measuring neutrino disappearance are of interest, as they constrain the ratio of the Dirac and Majorana 
masses. Additionally, we discuss electroweak constraints on C/(1)b-l, and show there is sufficient parameter space for 
a strong potential. 

We concentrate our analysis on fitting the positive anti-neutrino oscillation results of LSND and the null oscillation 
results of MiniBooNE above 475 MeV. We then discuss the speculative low energy excess seen at MiniBooNE and 
the ability of the model to accommodate this result. Finally, we give a prediction for anti-neutrino oscillations at 
MiniBooNE. 



V. CONSTRAINTS ON PARAMETERS FROM NEUTRINO OSCILLATION EXPERIMENT 

Before undertaking the numerical analysis, we consider constraints on parameters for the oscillation probability 
from experiment. The parameter space is large, and we make some simplifying assumptions. 

If the eigenvalues of the Dirac matrix are exactly equal, then in the small m/M limit, the matrix u will be the 
matrix governing mixing in the active sector as given by solar and atmospheric oscillation experiments in ([2|. We 
make this assumption, and pick the most straightforward values for ^13 and ^23 suggested by these experiments, 
013 = and 6*23 = it /A. We use the central value for the remaining angle, tan^ 9i2 ~ 0.45. 



The Dirac and Majorana neutrino masses m and Mi are loosely bound by experiments. The Dirac neutrino mass 
sets the scale for the Majorana neutrino masses, since their relative size governs the mixing between the active and 
sterile sectors. We take the Dirac mass to be sufficiently small, while still maintaining degeneracy in the neutrino 
masses, in the range m = 0.1 — 0.4 eV. Tests for active neutrino disappearance will constrain the Majorana masses 
Mi based on the range for m. In general, the Majorana mass will be in the range 1-10 eV. 

Electron neutrino disappearance tests give the most stringent constraint on Mi relative to m, as does the measure- 
ment of the total active solar neutrino flux at SNO |66j . Muon neutrino disappearance and conversion experiments 
[67l [68l l69l [70| FZll l72] involve baselines through matter at higher energies where the potential term suppresses the 
effect of the sterile neutrinos, so they will not provide a constraint. The best electron neutrino disappearance con- 
straint comes from the CHOOZ experiment [731 El]. Anti-electron neutrinos at CHOOZ have an energy of a few 
MeV, and propagate approximately 1 km through air before detection. The low B-L charge density in air coupled 
with the low neutrino energy means the potential term has essentially no effect. The SNO results are less sensitive to 
the disappearance of electron neutrinos than is the CHOOZ experiment. 

Both experiments are sensitive to oscillations with neutrino mass squared differences near Amg„„, far below the 
scales m and Mi relevant for oscillations between active and sterile neutrinos. Therefore, the oscillatory terms average 
out, and we can replace the oscillatory sin^ terms with i. The rate of electron neutrino disappearance into sterile 
neutrinos is therefore given by 

P{ve -^Vs)^ 2c?2 cos^ 6*1 sin^ 6*1 -I- 2si2 cos^ 62 sin^ 6*2 (24) 

where 



M? '''^\M^ 



sm^e.^^+of^) (25) 



SNO provides a measure of the total neutrino flux that agrees well with solar models, with a 9% experimental error 
and 16% theoretical error in the solar model [75] CHOOZ provides a better constraint in this case, having found 
no evidence for disappearance with a 3.9% total uncertainty. This gives a constraint 

9 ? 

TIT TTl 

0.039 > 1.37^-^0.62^ (26) 

The interference between oscillations involving the different eigenstates should average to zero, unless the masses are 
extremely degenerate. We have separate absolute bounds on Mi and M2: 

^>6 ^>4. (27) 

m m 

Since the sterile and active sectors decouple as Ali/m — s- 00, the only upper bound on this ratio will be set by LSND 
and MiniBooNE, when the sterile sector becomes relevant. 

The scale for the potential V is set by a number of considerations. Since it becomes relevant for LSND and 
MiniBooNE, in range of MiniBooNE energies the potential term 4VE should be comparable to or larger than Mf. 
Specifically, for oscillations at higher energies at MiniBooNE to be suppressed, we require the limit 

At LSND, however, oscillations should not be highly suppressed, at least not for anti-neutrinos. Therefore the potential 
should not be large compared to the quantity Mf / (AE^snu) ■ Since i^LSNo ~ O.I-Bmb, this implies 

<V< ' . (29) 

4i?MB ^ ^ 4_ElsnD 

Thus the optimal range for V is determined by the neutrino energy at the LSND and MiniBooNE experiments and 
the range of Majorana masses Mi. We cannot have V^^lsnd too large, and the range of Mf should be consistent with 
the LSND analysis and KARMEN constraint |77j, in the eV range. Therefore, we have bounds for the parameters 
that control the oscillation probability at LSND and MiniBooNE. After describing constraints on the vector mass and 
gauge coupling of C/(1)b-l! which constrain the potential V, we go on to describe our numerical analysis. 



Constraints on Weakly Coupled [/(1)b-l 



Precision electroweak measurements and rare decays constrain the coupling and vector mass for a U{1) gauge 
interaction [78l[79l|80l[81j[82j|83l[84l|85l|86]. Collider experiments can only constrain B - L if the couphng is 



relatively large, and the parameters we consider here are far away from this region. The most stringent constraint 
arises from corrections to g — 2 of the electron and muon. To leading order in the coupling, the electron g — 2 correction 
goes as 
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(30) 



where m„ is the paraphoton mass and g is the gauge coupling. The correction to the muon is the same, but with 



The experimental bounds on corrections to the electron and muon g — 2 are [871 1 
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1.5 • 10-^2 



Ja;^^ = 1.2 • 10~ 



(31) 



The potential neutrinos experience in matter is generated by low momentum transfers via the B — L photon between 
charges in matter and neutrinos. The strength of the potential is 



V = 



9 



rPB-L 



(32) 



where Pb l is the B — L charge density in matter, typically about 10'^'' charges/m'^ for rock at the Earth's surface. 
Since we want potentials stronger than 0.1 neV, we have a bound on .g/m„: 



— > 1.14- 10"^° eV"^ 



10- 



> 



88 keV 



(33) 



Since g must be small to meet the bounds from the electron and muon g ~ 2, we need my < m^. Then the best 
constraint comes from the electron g — 2, which, when combined with the independent measurement of a from the 
quantum Hall effect [H^, bounds g: 



g < 1.7-10" 



(34) 



To have potentials stronger than 0.1 neV, with g satisfying this bound, constrains m^,: my < 150 keV. If we use the 
lower limit of the potential found in numerical analysis, V > 3 neV, and thus the vector mass must be lighter than 
28 keV. 

For such a weak coupling, constraints on the new boson from rare particle decays are well within the bounds. The 
rate for any process with a paraphoton replacing a photon, such as e+e^ — s- 77*, is suppressed by a factor of g^/e^, 
which from the limit above, is 



^ <3-10-9 



(35) 



which is below the measured limits on rates of rare processes. 



VI. ASTROPHYSICS 



Note that the interactions of light gauge bosons with quarks, leptons and neutrinos, as well as the oscillations of 
sterile neutrinos, are severely constrained from their contributions to energy loss and energy transfer in stars and 
supernovae [53 We can avoid the stellar constraints if the effective paraphoton mass inside the star is much 
heavier than the temperature. For constant mass paraphotons these constraints pose a challenge to obtaining a 
sufficiently large potential to fit all the neutrino oscillation data. However the environment in the core of a star or 
supernova is extreme, with densities of B — L and neutrinos much higher than in the earth. In the next section we 
show that in a dense large object the effective mass of the paraphoton can be larger than it is on the earth, increasing 
the parameter space for paraphotons. 



A. Chameleon effect in stars 



Note that we have assumed that the paraphoton mass is due to the Higgs mechanism from a scalar s with B — L 
charge of 2. The field s is neutral under the standard model gauge group. The terms in the Lagrangian involving this 
scalar and the paraphoton field are 
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(36) 



Note that s has no direct couphng to ordinary matter or the active neutrinos. However ordinary matter has a net 
B — L charge density Pb l, which acts as a source for the Bq field. The Bq field in turn acts as negative term in 
the potential for s, increasing its expectation value and the effective mass of the paraphoton. Of course the Bq field 
is not gauge invariant, and neither is the phase of the s field. However the screening length, or inverse mass of the 
paraphoton, is a gauge invariant physical quantity which we wish to determine. 

To determine the paraphoton mass inside large objects, we consider a time independent spatially constant back- 
ground charge B — L charge density Pb-l- The s field carries a charge density ps which cancels the background 
charge. 



Ps 



-i2{Df)Ss* — sDqs*) 



(37) 



For a minimal energy configuration of given charge, we take s to have the form |s|e *2iot^ Since the phase of s is not 
gauge invariant it is convenient to define a gauge invariant field 



V = w + gBo 



(38) 



The charge density carried by the s condensate is — 4y|sp while V is the condensate energy per unit charge, which 
is identical to the neutrino potential energy used in the previous sections. 

To minimize the total energy for a static configuration where V and \s\ are spatially constant, one must solve the 
equations 



4:V\.S\ 



(39) 



The solution is 



m^2i -3 



V = 



16m° 
27 



- As/0 



2hUiXsPr + ^~XsP 



(40) 



and 



In the limit 



an approximate solution is 



|m| < (Ap)3 



(41) 



(42) 



V 



Note that the effective paraphoton mass is -v/S^jsl 



16 
P 

2^/x: 



(43) 



The chameleon effect is significant when the value of A satisfies Eq. ( 42 1 . The constraint in turn bounds the mass 
of the physical scalar excitations of the condensate. 
In the opposite limit 



|m| > (Asp)3 

the chameleon effect is negligible and an approximate solution is 

Am? 
\s\ w m/^/Xs 



(44) 



(45) 
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B. Numerical estimate of chameleonic paraphoton effect on stellar energy loss 



In this section we give a numerical example demonstrating how the chameleon effect can suppress the energy loss in 
stars to paraphotons. This is important because constraints from stellar evolution require that at densities between 
0.6 X 10'*g/cm'^ and temperatures around 0.7 x lO^K, and also for higher density and temperature of 2 x lO^g/cm^, 
and 10* K, radiation of exotic particles should not lead to energy loss greater than about 10 ergs/(g-s)[55]. 

For simplicity in the analysis of stellar constraints we take g = 10~^, niy = 28 keV. In order to have a significant 
chameleon effect at stellar densities we take a very flat potential for the scalar field, rus = 10~® eV, and A<; = 5 x 10~^^. 
(Note that because ordinary matter does not act as a source for the scalar field, such a light s field is consistent with 
constraints from searches for new forces.) The paraphoton mass at typical earth densities of 2g/cw? is then nearly 
the same as in vacuum, while at a density of 0.6 x 10*g/cm^, assuming about 1/2 of the baryons are neutrons, the 
paraphoton mass is 260 keV. A typical temperature for stellar cores of such a density is 0.7 x 10* K and so the 
Boltzmann suppression factor for paraphoton emission is about 10~^^. The emission rate for paraphotons would then 
be similar to the rate for emission of light vector bosons with an a of order 10^"^", which is smaller than the quoted 
stellar evolution limit 'SS] for that density of 10~^^. At a density of 2 x lO^g/cm'' the paraphoton mass is 800 keV, 
and the Boltzmann suppression factor at a temperature of 10* K is completely negligible. 



C. Supernova 1987A 



Weak interactions of bosons lighter than ~ 100 MeV are constrained by the observation of a small number of 
neutrino events lasting for about 10 seconds after the explosion of SN 1987A |.89. .90. 91 . If the paraphoton is much 
lighter than the supernova temperature, couplings stronger than about 10~^ are excluded, as neutrinos would then 
remain thermally coupled until lower temperatures than suggested by the energies detected [52], while couplings 
weaker than about 10~^ are also excluded, as then the supernova would cool too rapidly by emitting paraphotons 

[Miins]- 

For the parameters considered in SVIB at a typical nuclear density of lO^^g/cm'^ found in a supernova core, the 
paraphoton mass would be of order 700 GeV. Even if the potential is not this flat, and so the paraphoton is lighter 
than the temperature and in thermal equilibrium with the neutrinos, it is strongly enough coupled to not lead to 
excessively rapid supernova cooling . 

One also must consider the question of neutrino emission. Usually the supernova observation provides severe 
constraints on sterile neutrinos because they are more weakly coupled and can more easily escape than neutrinos. 
The sterile neutrino interaction cross section depends on the paraphoton mass. If the paraphoton is light the sterile 
neutrino cross section from paraphoton exchange at a temperature below 100 MeV is larger than the weak cross 
section, and so sterile neutrinos would not escape more easily than ordinary neutrinos. For a paraphoton mass of 700 
GeV the sterile neutrino productions cross section is of order 10~^* times the weak cross section and so production 
of sterile neutrinos may be neglected. A range of intermediate high density paraphoton masses will be excluded by 
supernova cooling, however astrophysical constraints are not the main focus of this paper and so this range will not 
be computed here. 

One also must consider sterile neutrino production via neutrino oscillations. However, for such large values of the 
B-L potential as found in supernova core, active sterile neutrino mixing angles are negligible. 



D. Nucleosynthesis constraints 



The mechanism used for baryogenesis will determine the net B — L charge of the universe, and the B — L charge 
density in the early universe. Electroweak baryogenesis, for instance, gives a net B — L charge of zero while 
Affleck-Dine baryogenesis [95] could give a large net B ~ L density. 

If we assume a net B — L charge of zero, then the main change to the effective potential for s comes from thermal 
corrections, including the neutrino mass contribution to the energy density. Even with a potential as flat as assumed 



in section SVIB these only increase the s vev by a factor of about a 100, and so the paraphoton and the sterile 
neutrinos are lighter than the temperature during nucleosynthesis. The paraphoton may be produced singly and will 
be in thermal equilibrium, acting approximately as another neutrino species. The cross section for sterile neutrino 
interactions at nucleosynthesis temperatures from paraphoton is comparable to weak, and so the total number of 
effective neutrino species during nucleosynthesis would be 7, a value inconsistent with the Helium abundance j89j . 

Another attractive mechanism, available in supersymmetric theories, is Affleck-Dine baryogenesis, which relies on 
large scalar field vevs in the early universe along flat directions of the scalar potential. This mechanism could create 
a large net B — L density along a flat direction where s has a large vev. The minimum of the effective potential 
for s would then decrease with the Hubble expansion. This has three potential effects on nucleosynthesis. First, 
a large lepton number chemical potential would affect the Helium abundance, and allow for more light relativistic 
species [SSI IHT] Secondly, even it the lepton number chemical potential is not large compared with the nucleosynthesis 
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temperature, the paraphoton and sterile neutrino masses could always be larger than the temperature, and so these 
are never in equilibrium. Third, even if the sterile neutrinos are lighter than the temperature, a heavier paraphoton 
reduces the cross section for sterile neutrino production, and so they need not be in thermal equilibrium. Note that 
sterile neutrino production via neutrino oscillations will be negligible due to tiny active-sterile mixing angles at high 
B — L charge density and large s vev. 



VII. NUMERICAL ANALYSIS 



We do a numerical analysis of the oscillation probability at LSND and MiniBooNE to find the region of parameter 
space that fits both of these experiments. We scan over parameters that satisfy the CHOOZ electron disappearance 
bounds. For these parameters the bounds on muon neutrino short baseline disappearance |S7] are also satisfied. We 
then make a quantitative prediction for the results of MiniBooNE anti-neutrino oscillations. We are interested in 
fitting the anti-neutrino oscillation results of LSND and the null oscillation results above low energy at MiniBooNE. 
The low energy excess seen at MiniBooNE and anti-neutrino oscillations at MiniBooNE may then be considered in 
the context of these two results. The results can be summarized in terms of the oscillation probability: 

PlsndC^^m ^ ^e) = (0.264 ± 0.081)% 
-PmbIj^m ^ 475 MeV <E <i GeV) < 0.064% 

For each experiment, the energy distribution of neutrinos must be taken into account. We use the Michel spectrum 
from the muon decay giving rise to at LSND, and fit an interpolating function to the approximate energy 
spectrum at MiniBooNE [98 . If s{E) is the relevant energy density function, then the oscillation probability at a 
baseline L is just 

P{L) = j dEs{E)P{L,E) 

and we do not consider detector geometry at either experiment. 

We pick the following values for the parameters, and study the regions listed below for the remaining ones: 

6*12 0.591 = 023 = 7r/4 

0.1 neV < y < 5 neV 0.1 < to < 0.4 eV 
6to < Afi < 7 eV 4m < Ma < 7 eV 

9 9 

m TO 
0.039 > 1.37— ^-f 0.62— . 

The last limit above imposes the CHOOZ constraint in ref. [T^ . We choose an upper bound on V to place a reasonable 
lower bound (10 keV) on the paraphoton mass to„. The lower limit on V and the upper limit on cut off the analysis 
region where changes in the magnitude of those variables have a small effect on the oscillation probability. As we will 
see, those limits do not serve as real constraints on the parameter space, with very few points saturating the limits. 



(46) 



(47) 



(48) 



Results 



For the analysis, we sample the parameter space with 3 • 10^ points. At each point, we require the LSND anti- 
neutrino result and the MiniBooNE null result to be satisfied, as well as the CHOOZ constraint. For the points that 
pass the constraints, we find both the MiniBooNE low energy oscillation probability and the MiniBooNE anti-neutrino 
oscillation probabilities at high and low energy. The total volume of the parameter space we fit is 9.95 eV'^- 4.9 neV, 
giving a mean volume spacing between points of (0.05 eV)^- 0.11 neV. Because the oscillation probability is smooth, 
we are confident that this coverage is sufficient to find interesting regions in the parameter space. 

Approximately 0.1% of the trial points fit LSND and MiniBooNE and pass the CHOOZ constraint, meaning there 
is an ample region of parameter space which reconcile LSND and MiniBooNE. Of the points fitting the constraints, 
the parameters fall within the ranges: 

2.5 neV <V <5 neV 0.2 eV < m < 0.35 eV 
2.5eV<Mi<7eV 0.9 eV < M2 < 1.7 eV 

Many points sit near the CHOOZ bound, implying it provides a real constraint on the model. Conversely, the lower 
limits for m and V, and the upper limit on the Majorana masses M^ that we have chosen do not seriously constrain the 
ability of the model to reconcile LSND and MiniBooNE, as few points in the fits lie close to these bounds. Therefore, 
we conclude that the parameter space we have chosen is suitable to test whether the model can reconcile LSND and 
MiniBooNE. 



(49) 
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The Low Energy Excess at MiniBooNE 



For the points passing the fits, we find the model is capable of fitting a portion of the observed low energy excess at 
MiniBooNE. If the low energy excess observed at MiniBooNE is determined to be a signal of — > transmutation, 
then we are obviously interested in fitting such a dramatic excess with this model. However, we find difficulty fitting 
the entire excess while maintaining the observed low values of the oscillation probability at LSND and MiniBooNE 
(above 475 MeV). In the energy range of 200 to 475 MeV, we fit at most approximately 40% of the observed excess, 
corresponding to an oscillation probability of 0.5%. A number of points have a small low energy excess, which are 
suitable to reconcile LSND and MiniBooNE if the observed MiniBooNE excess reduces with further analysis. 

We have restricted the parameter space (see eq. 48 1 in order to make the numerical analysis faster. In particular 



we do not vary the mixing angles in the matrix u, which control the flavor fractions of the active component of the 
heavier neutrinos. In our analysis, only 2 of the heavy neutrinos have any electron flavor component. One effect of 
varying the angles in u would be to allow all three sterile neutrinos to play a role. Allowing one of the heavy neutrinos 
to have a small electron flavor component would allow it play a role in electron flavor conversion while being lighter 
than 4m without violating the CHOOZ constraint. Should the low energy MiniBooNE electron excess survive further 
analysis, it would be interesting to see if the model can give the entire excess when more parameters are varied. 
However we believe we have varied enough parameters to cover the essential implications of the model. 

In Figure 1, we give a histogram of the low energy excess for the flt points, and in Figure 2 we show the correlations 
between the low energy excess and the oscillation probability above 475 MeV at MiniBooNE. Figure 2 shows the 
main difficulty in fitting the low energy excess in the model — there is a correlation between the well constrained 
MiniBooNE oscillation probability above 475 MeV and the low energy excess. 



MiniBooNE Oscillation Probability at Low Energy 




FIG. 1: A histogram of tlie MiniBooNE oscillation probability for neutrinos in the energy region from 200 to 475 MeV 



0.005 - 



0.002 - 



MiniBooNE Oscillation Probability 




0.0002 0.0004 

Above Low Energy 



FIG. 2: The average MiniBooNE neutrino oscillation probability in the energy region from 200- 475 MeV vs. the region above 
475 MeV 
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Predictions for Anti-Neutrino Oscillations at MiniBooNE 

We can make definite predictions for anti-neutrino oscillations at MiniBooNE. Because |yi?MB| must be similar to 
the Majorana masses Mf to reconcile LSND and MiniBooNE, in anti-neutrinos we expect an MSW type resonance 
near or in the low energy MiniBooNE region which leads to relatively large — > v^. conversion. At our fit points, 
we can examine the oscillation probability in ant i- neutrinos and correlate the predicted signal with the oscillation 
probability in neutrinos at MiniBooNE. In our analysis, we assume that the energy spectrum of anti-neutrinos is the 
same as neutrinos at MiniBooNE. In Figure 3, we give the anti-neutrino oscillation probability in the low energy 
region from 200- 475 MeV versus the region above 475 MeV. 



MiniBooNE Anti-Neutrino Oscillation Probability 

0.25 r 




0.02 0.04 0.06 0.08 0.1 

Above Low Energy 

FIG. 3: The MiniBooNE anti-neutrino oscillation probability in the region from 200 to 475 MeV vs. the region above 475 MeV 



Due to the resonance in anti-neutrinos in or just below the MiniBooNE energy range, the anti-neutrino oscillation 
probability is consistently much larger than the neutrino oscillation probability over the energy range of the experiment. 
We can express this quantitatively over the fit points by the ratio between the anti-neutrino and neutrino oscillation 
probabilities. In each energy region, the minimum ratio is above 3; in Figure 4, we give the cumulative distribution 
of the fraction of fit points with an anti-neutrino to neutrino oscillation probability ratio above a fixed value. 



Cumulative MiniBooNE Anti-Neutrino to Neutrino Ratio 




10 20 30 40 50 

Oscillation Probability Ratio 

FIG. 4: The cumulative distribution for ratios of anti-neutrino to neutrino oscillation probability, below (black, upper curve) 
and above (blue or gray, lower curve) 475 MeV energy. For a fixed oscillation probability ratio R, the vertical axis is the fraction 
of fit points with ratio greater than R 



Our analysis leads us to make a quantitative prediction for MiniBooNE: anti-neutrino oscillations are likely to be 
observable at the experiment. The experimental systematic uncertainty may become relevant for for fit points that 
predict particularly low values of the anti-neutrino oscillation probability in both energy regions, but the majority 
of fit points correspond to anti-neutrino oscillation probabilities that are large enough to be free of this concern, as 
shown in Figure 3. 
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VIII. SUMMARY 



We have shown at a model with three sterile neutrinos and a light _B — L gauge boson can accommodate the results 
of both LSND and MiniBooNE. which are in direct conflict under a standard oscillation model. The effects of the 
light vector boson with weak gauge coupling can be made to fit experimental limits, while having a significant effect 
on short baseline neutrino oscillations. We have shown agreement with other neutrino oscillation experiments besides 
LSND and MiniBooNE, and made a quantitative prediction that anti-neutrino oscillations at MiniBooNE are much 
larger than neutrino oscillations, and, with sufficient anti-neutrino data, are likely to be visible. 
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